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1 Random Matrix Models 

In probability theory and statistics a common first approximation to many random processes 
is a sequence Xi, X2, X3, ... of independent and identically distributed (iid) random variables. 
Let F denote their common distribution. To motivate the material below, we take these random 
variables and construct a particularly simple N x N random matrix, 

diag {Xi{uj), X2{uj), ■ ■ ■ , Xn{uj)) ■ 

The order statistics are the eigenvalues ordered 

Ai < A2 < • • • < Aiv, 

and the distribution of the largest eigenvalue, Amax(-^) = Aat, is 



Since the distribution F is arbitrary, we observe that so too is the distribution of the largest 
eigenvalue of a A'^ x random matrix. However, one is really interested in limiting laws as 
— > 00. That is, we ask if there exist constants a^v and such that 



Prob(Amax(iV) <x) 



Prob {Xi<x,...,Xn <x) 
F{xf. 



converges in distribution to a nontrivial limiting distribution function G. In the present situa- 
tion a complete answer is provided by 

Theorem: If ( |1.1| ) converges in distribution to some nontrivial distribution function G, then G 
belongs to one of the following forms: 

1. e"*^ with support R. 

2. e'^^^"' with support [0, oo) and a > 0. 

3. e^("^)° with support (— oo,0] and a > 0. 

This theorem is a model for the type of results we want for nondiagonal random matrices. 

A random matrix model is a probability space {Vl,V,J-) where is a set of matrices. Here are 
some examples 

• Circular Unitary Ensemble (CUE, (3 = 2) 

— Vl = U{N) = N X N unitary matrices. 

— 'P= Haar measure. 

• Gaussian Orthogonal Ensemble (GOE, /3 = 1) 

— Q = N X N real symmetric matrices. 

— V = uniqueQ measure that is invariant under orthogonal transformations and the 
matrix elements (say on and above the diagonal) are iid random variables. 

• Gaussian Unitary Ensemble (GUE, /3 = 2) 

— Q = N X N hermitian matrices. 

— V= unique measure that is invariant under unitary transformations and the real and 
imaginary matrix elements (say on and above the diagonal) are iid random variables. 

• Gaussian Symplectic Ensemble (GSE, (3 = 4) 

— Q = 2N X 2N Hermitian self-dual matrices .0 

— V= unique measure that is invariant under symplectic transformations and the real 
and imaginary matrix elements (say on and above the diagonal) are iid random 
variables. 

^Uniqueness is up to centering and a normalization of the variance. 

^Identify the 2 A*' x 2N matrix with the N x N matrix whose entries are quaternions. If the quaternion 
matrix elements satisfy Mji = Mij where the bar is quaternion conjugation, then the 2N x 2N matrix is called 
Hermitian self-dual. Each eigenvalue of a Hermitian self-dual matrix has multiplicity two. 
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Expected values of random variables f : ^ C are computed from the usual formula 

En{f)= f f{M)dV{M). 
If f{M) depends only on the eigenvalues of M G then one can be more explicit: 
• CUE (Weyl's Formula) 



EumU) = r ■■■T m, ...,on) n ' • • • ' 



2 

dOi ■ ■ ■ dO^i 



• Gaussian Ensembles {13 = 1,2,4): 



ENfsif) = CNf3 I ■ I f{xi,...,XN)\Mxi,...,XN)f e ^ T, dxi ■ ■ ■ dx N , 



where cn/3 is chosen so that E]y/^{1) = 1 and A{xi, . . . , xn) = Yli<i<:j<Ni^i ~ ^j)- "^^^ 

factor e 2 explains the choice of the word "gaussian" in the names of these ensembles. A 
commonly studied generalization of these gaussian measures is to replace the sum of quadratic 
terms appearing in the exponential with '^V{xi) where V is, say, a polynomial (with the 
obvious restrictions to make the measure well-defined). 

Choosing / = Hi (l ~ Xji^-iSj: XJ the characteristic function of a set J C R, we get the 
important quantity]^ 

En/sH) = E]\fi3{0; J) := probability no eigenvalues lie in J, 

and in the particular case J = {t,oo) 

FNp{t) := Prob (Amax < t) = E^piO, J). 

The level spacing distribution^ is expressible in terms of the mixed second partial derivative of 
En (3(0] (fl) b)) with respect to the endpoints a and b. 

2 Fredholm Determinant Representations 

Though E]\fi3{0;J) are explicit A^-dimensional integrals, these expressions are not so useful 
in establishing limiting laws as — > 00. What turned out to be very useful are Fredholm 

^This quantity has an obvious extension to other random matrix models. 

*Let the eigenvalues be ordered. The conditional probability that given an eigenvalue at a, the next one lies 
between s and s + ds is called the level spacing density. 
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determinant representations for Ei^p{{)]J). In 1961 M. Gaudin proved for [5 = 2 (using the 
newly developed orthogonal polynomial method of M. L. Mehta) that £^Ar2(0; J) = det(/— K7V2) 
where Kj^2 is an integral operator acting on J whose kernel is of the form 

^{x)ip{y) -i'{x)ip{y) 
x-y ' 

with ip{x) = C]\[e~^^^^H]\[{x), il}{x) = CNe~^^^^H]\f-i{x), and Hj{x) are the Hermite polyno- 
mials.0 For /3 = 1 or 4, generalizing F. J. Dyson's 1970 analysis of the n-point correlations for 
the circular ensembles, it follows from work by Mehta the following year that the square of 
EN/si^; J) again equals a Fredholm determinant, det(/ — K^p)^ but now the kernel of K^p is 
a 2 X 2 matrix.0 



3 Scaling Limits (Limiting Laws) 
3.1 Bulk Scaling Limit 

Let pn{x) denote the density of eigenvalues at x and pick a point xq^ independent of A'" with 
Pn{xo) > 0. We scale distances so that resulting density is one at xq, ^ := pwixo) {x — xq), 
and we call the limit 

— > 00, x — > Xq, such that ^ is fixed, 
the bulk scaling limit. For (3 = 2, 

En2{0; J) ^2(0; J) = det (/ - K2) 
where the integral operator K2 (acting on L'^{J)) has as its kernel (the sine kernel) 

1 sin7r(.^ — 

(We use the same symbol J to denote the scaled set J.) Furthermore, 

P2{s) = -^E2iO; (0,s)) 

is the (limiting) level-spacing density for GUE; known as the Gaudin distribution .[^ We observe 
that the limiting kernel is translationally invariant and independent of xq. 

^For the random matrix models corresponding to general potential V, ip{x) = Cffe~^^^^^^pM{x) and i/'(a;) — 
CNe~^^^^^^PN-i{x) where Pj{x) are the orthogonal polynomials associated with weight function w{x) = e~^'^'. 
It is in this generalization that we see the close relation between the general theory of orthogonal polynomials 
and random matrix theory. 

^ See for elementary proofs of these facts. 

'^For the analogous P — 1,4: results, see, e.g., or [Q. 
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3.2 Edge Scaling Limit 

In the gaussian ensembles, the density decays exponentially fast around 2a^/N ; perhaps sur- 
prisingly, it is also the case that 

Amax(^) „ /„ ,x 

lim = — = 2a, a.s. o.l) 



where a is the standard deviation of the off-diagonal matrix elements. (In the normalization 
we've adopted, a = l/\/2.) If we introduce the scaled random variable A through 

Amax = 2fT\/iV + 



iVl/6 ' 
then 

Prob (Amax < t) = Prob (a < s) ^ Fp{s) as TV ^ oo, 
where t = 2a Vn + as/N^/^. For /3 = 2, 

F2(s) = det(/-/CAiry)' 
where -f^Airv kernel of the form ( |2.1[ ) with (p{x) = Ai(x), ^(x) = Ai'(x) and J = {s, oo) 
(See, e.g., p% for the /? = 1,4 results.) 



4 Connections with Integrable Systems 
4.1 Bulk Scaling Limit 

In 1980 M. Jimbo, T. Miwa, Y. Mori, and M. Sato |21] expressed the Fredholm determinant of 
the sine kernel in terms of a solution to a certain system of integrable differential equations.^ 
In the simplest case of a single interval, J = (0, s), the differential equation is a particular case 
of Painleve V (-Py)^ and the Fredholm determinant is given by 

det (/ -XK2) = exp (^J^' dx^ , 

{xa")^ + 4 (xa' - a) (xa' -a + {a'f^ = 0, 

. XX A 
a{x; A) ~ X, as x — > 0. 

TT 

For /3 = 1,4 and J = (0, s), Ef^{0; (0,s)) can also be expressed in terms of the same function 
a{x; 1). A down-to-earth application of these Painleve representations (and using the known 
asymtptotics of a{x]l)) is that one can easily produce graphs of the level spacing densities 

ppjsm 

*A simplified proof of tlieir results can be found in [Q. 

''The differential equation below is the sigma representation of Pv- 

^"Without the Painleve representations, the numerical evaluation of the Fredholm determinants is quite in- 
volved. 
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Figure 1: Densities for the scaled largest eigenvalues, fpis). 
4.2 Edge Scaling Limit 

The limiting distributions (edge scaling) of the largest eigenvalue, Ff^{s), can also be expressed 
in terms of Painleve functions — this time Pjj [3£, 40 1: 



Fi(s)2 = exp 
F2{s) = exp 
F4{s/V2f = cosh 
where q satisfies the Painleve II equation 



q{x)dxj F2{s), 
{x — s)q{x)'^ dx^ , 
q{x)dx] F2{s), 



(4.1) 
(4.2) 
(4.3) 



q" =xq + 2q^ 



with boundary condition q{x) ~ Ai(x) as x ^ oo.0 The graphs of the densities //3(s) 
dFp{s)/ds are in Figure 1. 



^^That a solution q exists and is unique follows from the representation of the Fredholm determinant in 
terms of it. Independent proofs of this, as well as the asymptotics as x — cxd, were given by S. Hastings and 
J. McLeod, P. Clarkson and McLeod and by P. Deift and X. Zhou. 
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4.3 Generalizations 



Both the sine kernel and the Airy kernel are of the form (|2.1j ). Kernels of this form arise in 
many problems in integrable systems; indeed, so much so that A. Its, A. Izergin, V. Korepin and 



V. Slavnov in 1990 initiated a general analysis of these kernels. The following theorem [41|, 
which applies to a wide class of /? = 2 random matrix ensembles, gives the general situation: 

Theorem: Let J = UjLi('^2j-i — CL2j) be a union of open intervals. Define T(a) = det (/ — K) 
where K is an integral operator acting on Lp'{J) whose kernel is of the form ( |2.1D where ip and 
tp are assumed to satisfy 

with a 2 X 2 matrix with zero trace and rational entries in x, then logdet(I — K) are 
expressible polynomially in terms of solutions to a total system of partial differential equations 
{ttj are the independent variables). The differential equations are given explicitly in terms of 
the coefficients of the rational functions appearing in n{x). 



4.4 Historical Comments 

The first connection between Toeplitz/Predholm determinants and Painleve functions was es- 
tablished in 1973-77 in work of T. T. Wu, B. M. McCoy, E. Barouch and the first author 
concerning the scaling limit of the 2-point functions of the 2D Ising model of statistical me- 
chanics. The Painleve function that arose was Puj. This work was subsequently generalized 
by Sato, Miwa and Jimbo to n-point functions and, more generally, holonomic quantum fields. 
The Kyoto School then took up the problem of the density matrix of the impenetrable Bose 
gas and it was in this context that they discovered that the Fredholm determinant of the sine 
kernel is related to Py- 

A crucial simplification of the Kyoto School work, as it applies to random matrix theory, was 
made by Mehta in 1992 [^8|. This last work inspired the commutator methods introduced by the 
present authors in the period 1993-96. Since then both Riemann-Hilbert methods of Deift, Its, 
Zhou and others (see, e.g. [|l5|); and Virasoro methods of M. Adler, T. Shiota, P. van Moerbeke, 
and others (see, e.g. [||), have played an increasingly important role in the development of 
random matrix theory. The connection of these methods with the isomonodromy method has 
been clarified by J. Palmer |^2| and J. Harnad |19|. 



Space does not permit us to discuss the interesting connections between random matrices and 
Szego type limit theorems. See E. Basor [R] for connections with linear statistics and the review 



papers 42 1 for some related historical comments. 
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5 Universality 



5.1 Universality of Gaussian Ensembles in Random Matrix Models 
5.1.1 Invariant Measures, (3 = 2 

As briefly mentioned above, a widely studied class of random matrix models is defined by 
the replacement of the gaussian potential, , by the general potential V{x). For the weight 
functions most studied, the parameter N is put into the exponent so that the weight function 
becomes ^-^'^i^) _ For different F's, the limiting density pv{^) can be quite different. It may be 
supported on many distinct intervals, and it may vanish at interior points of its support. In the 
gaussian case, the limiting density is the Wigner semicircle law: pw{x) = :|Vl — x"^. Heuristic 
arguments suggest that the behavior exhibited by the Wigner law — that p is positive on the 
interior of its support and vanishes like a square root at endpoints — is the typical behavior 
for py. The bulk scaling limit and edge scaling limit are defined in analogous ways to the 
gaussian cases. To establish universality of these scaling limits, one must show (for (3 = 2 
ensembles) that the scaled kernels approach the sine kernel and the Airy kernel, respectively. 
The potential V{x) = ^x^ + fx^ {g > 0, t < 0) is an example of a "two interval" potential. 
Indeed, for this important example P. Bleher and A. Its Q proved precisely this statement of 
universality. (See their paper for related work in the orthogonal polynomial literature as well as 
the physics literature.) Recently, building on work of fl^, A. Kuijlaars and K. McLaughlin |26| 



have shown this behavior is generic for recti analytic Scitisfying lini^ y(x)/log|x| = +00. 

In the physics literature, M. Bowick, E. Brezin [|T2| and others have argued (for (5 = 2 ensembles) 
that if pv vanishes faster than a square root, then the corresponding edge scaling limit will 
result in nonAiry universality classes. The resulting new kernels will have form (|2.1| ) and the 
theory developed in [^] will apply, but there remains much to be understood in these cases. 



For /3 = 1,4, the situation is more complicated due to the structure of Ky |37|, and the 
"universality" theorems are not so general. 



5.2 Noninvariant Measures: Wigner Ensemble 

The Wigner ensembles are defined by requiring that the matrix elements on and above the 
diagonal in either the real symmetric case or the complex hermitian case are independent and 
identically distributed random variables. It is only in the case when the distribution is gaussian 
is the measure invariant. One usually assumes, as we do here, that all moments of the common 
distribution function exist. It was Wigner himself who showed that the limiting density of 
states is the Wigner semicircle. Subsequently several authors — culminating in a theorem by 
Z. Bai and Y. Lin clarifying which moments need exist — showed (^]^) continues to hold for the 
Wigner ensembles. 

It should be noted that because the measure is noninvariant, the nongaussian Wigner ensembles 
do not, as far as we understand, have Fredholm determinant representation for their distribu- 
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tion functions. This means, for one, that the methods of integrable systems are not directly 
apphcable to Wigner ensembles. It is therefore particularly important, as A. Soshnikov |36| 
recently proved, that in the edge scaling limit the Wigner ensembles are in the same universal- 
ity class as the gaussian models. In particular, the limiting distribution of the scaled largest 
eigenvalue is given by -Fi(s) for real symmetric Wigner matrices and by F2{s) for complex 
hermitian Wigner matrices. 



5.3 Examples from Physics 

A second type of universality, and the one first envisioned by Wigner in the context of nuclear 
physics, asserts in Wigner's words [46| 

Let me say only one more word. It is very likely that the curve in Figure I [an 
approximate graph of Pi{s)] is a universal function. In other words, it doesn't 
depend on the details of the model with which you are working. 

The modern version of this asserts that for a classical, "fully" chaotic Hamiltonian the corre- 
sponding quantum system has a level spacing distribution equal to pp{s) in the bulk. (The 
symmetry class determines which ensemble.) This quantum chaos conjecture, due to O. Bohi- 
gas, M. Giannoni and C. Schmit |1C], has been a guiding principle for much subsequent work, 



though it is the authors' understanding that it remains a conjecture. A particularly nice nu- 
merical example supporting this conjecture is M. Robnik's work [^] on chaotic billards. The 
reader is referred to the recent review article [^] for further numerical examples that support 
this conjecture. It should be noted that there are examples from number theory where the 
conjecture fails. Thus, as it has been said, the conjecture is undoubtedly true except where it 
is demonstratively false. 



5.3.1 Aperiodic Tiling Adjacency Matrix 

The discovery of quasicrystals has made the study of statistical mechanical models whose un- 
derlying lattice is quasiperiodic of considerable interest to physicists. In particular, in order to 
understand transport properties, tight binding models have been defined on various quasiperi- 
odic lattices. One such study by Zhong et al. [^] defined a simplified tight binding model 
for the octagonal tiling of Ammann and Beenker. This quasiperiodic tiling consists of squares 
and rhombi with all edges of equal lengths (see Figure 2) and has a symmetry around the 
central vertex. On this tiling the authors take as their Hamiltonian the adjacency matrix for 
the graph with free boundary conditions. The largest lattice they consider has 157,369 vertices. 
The matrix splits into ten blocks according to the irreducible representations of the dihedral 
group D^. For each of these ten independent subspectra, they compare the empirical distri- 
bution of the normalized spacings of the consecutive eigenvalues with the GOE level spacing 
density pi{s). In Figure 2 we have reproduced a portion of their data for one such subspectrum 
together with pi . 
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Figure 2: Data for nearest neighbor normalized spacings of eigenvalues of the adjacency matrix 
for a quasiperiodic octagonal tiling are plotted together with the GOE level spacing density 
Pi{s). Data are from one independent subspectrum of a Dg-symmetric octagonal patch of a 
tiling with 157,369 vertices. Courtesy of Zhong et al. [p7| ]. 



5.4 Spacings of the Consecutive Zeros of Zeta Functions 

Perhaps the most surprising appearance of the distributions of random matrix theory is in 
number theory. Analytical work by H. Montgomery and extensive numerical calculations by 
A. Odlyzko on the zeros of the Riemann zeta function have given convincing evidence that 
the normalized consecutive spacings follow the Gaudin distribution, see Figure 3. Recent 
results of Z. Rudnick and P. Sarnak are also compatible with the belief that the distribution 
of the spacings between zeros, not only of the Riemann zeta function, but also of quite general 
automorphic L-functions over Q, are all given by this Montgomery-Odlyzko Law. In their 
landmark book |25], N. Katz and P. Sarnak establish the Montgomery-Odlyzko Law for wide 
classes of zeta and L-functions over finite fields. 



5.5 Random Matrix Theory and Combinatorics 

The last decade has seen a flurry of activity centering around connections between combinatorial 
probability of the Robinson-Schensted-Knuth (RSK) type on the one hand and random matrices 
and integrable systems on the other. From the point of view of probability theory, the quite 
surprising feature of these developments is that the methods came from Toeplitz determinants, 
integrable differential equations of the Painleve type and the closely related Riemann-Hilbert 
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1 ,041 ,600 zeros near the 2*10'^20-th zero 
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Figure 3: Data for nearest neighbor spacings among 1,041,600 zeros of the Riemann zeta 
function near the 2 x 10^*'-th zero are plotted together with the GUE spacing density. Courtesy 



of Andrew Odlyzko 



techniques as they were apphed and refined in random matrix theory. Using these techniques 
new, and apparently quite universal, limiting laws have been discovered. The earliest signs of 



these connections can be found in the work of A. Regev |34| and I. Gessel |17]. Here, however, 
we introduce this subject by examining a certain card game of D. Aldous and P. Diaconis Q, 
called patience sorting. 



5.5.1 Patience Sorting and Random Permutations 

Our deck of cards is labeled {1, 2, . . . , N} and we order the cards with their natural ordering. 
Shuffle the deck of cards and 



Turn over the first card. 

Turn over the second card. If it is of higher rank, start a new pile to the right of the first 
card. Otherwise place the second card on top of the first card. 

Turn over the third card. If it is of higher rank than either the first or the second card, 
start a new pile to the right of the second card. Otherwise place the third card on top of 
the card of higher rank. If both first and second are of higher rank, place the third card 
on the smaller ranked card. (That is, play cards as far as possible to the left.) 
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Continue playing the game, playing cards as far left as possible, until all the cards are 
turned over. 



The object of the game is to end with a small number of piles. Let ^Ar(cT) equal the number 
of piles at the end of the game where we started with deck a = {ii,i2, ■ ■ ■ jIn}- Clearly, 
1 < iN{cT) < N, but what are some typical values for a shuffled deck? Starting each time with 
a newly shuffled deck of = 52 cards, the computer played patience sorting 100,000 times. 
Here are the statistics for £52: 



• Mean=11.56 (11.00). 

• Standard Deviation=1.37 (1-74) 

• Skewness=0.33 (0.22) 

• Kurtosis Excess =0.16 (0.09) 

• Sample Range = 7 to 19 (Probability 0.993) 



where the numbers in parentheses are the asymptotic predictions (as the number of cards tends 
to infinity) of the theory of J. Baik, P. Deift and K. Johansson Q to be described below. 

A shuffled deck of cards, a = {ii, ^2, • • • > ^Af}i is a permutation of {1, 2, . . . , n}, and so we think 
of the shuffled deck as a random permutation. A moment's reflection will convince the reader 
that In{cf) is equal to the length of the longest increasing subsequence in the permutation a. 
As a problem in random permutations, determining the asymptotics oi E{In) as N 00 is 
called Ulam's Problem. In the 1970's A. Vershik and S. Kerov and independently B. Logan and 
L. Shepp showed E{£]\f ) ~ with important earlier work by J. Hammersley. Hammersley's 
analysis introduced a certain interacting particle system interpretation. This was developed 
by Aldous and Diaconis Q who in 1995 gave a "soft" proof of this result using hydrodynamic 
scaling arguments from interacting particle theory. 

Introducing the exponential generating function 

Y^ProhiiN<n) — , 

N>0 

Gessel showed that it is equal to Dn{t), the determinant of the n x n Toeplitz determinant 
with symbol e^^^"*"^ (Recall that the i,j entry of a Toeplitz matrix equals the i — j Fourier 
coefflcient of its symbol.) It is in this work of Gessel and subsequent work of Odlyzko et 



al. 1 30 1 and E. Rains |3^, that the methods of random matrix theory first appear in RSK type 



problems .p^ 



Gessel does not mention random matrices, but in light of well-known formulas in random matrix theory 
relating Toeplitz determinants to expectations over the unitary group, we believe it is fair to say that the 
connection with random matrix theory begins with this discovery. See, however, Regev p4|. 
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Starting with this Toephtz determinant representation, Baik, Deift and Johansson [0], using 



the steepest descent method for Riemann-Hilbert problems [14|, derived a dehcate asymptotic 
formula for -D„(t) which we now describe. Introduce another parameter s and suppose that n 
and t are related by n = [2t^/^ + st^^^]. Then as t ^ oo with s fixed one has 



lim e-'Dn{t) = F2{s) 



where F2{s) is the distribution function ( |4.2| ). Using a dePoissonization lemma due to Johans- 
son |2^, these asymptotics led Baik, Deift and Johansson to the limiting law 



Since the work of Baik, Deift and Johansson, several groups have extended this connection 
between RSK type combinatorics and the distribution functions of random matrix theory. The 
aforementioned result is equivalent to the determination of the limiting distribution of the 
number of boxes in the first row in the RSK correspondence a <-> {P,Q). In the same 
authors show that the limiting distribution of the number of boxes in the second row is (when 
centered and normalized) distributed as the second largest scaled eigenvalue in GUE [^. They 
then conjectured that this correspondence extends to all rows. This conjecture was recently 
proved by A. Okounkov [p^] using topological methods and by A. Borodin, A. Okounkov and 



G. Olshanski |13] and Johansson [^] using analytical methods. 

Placing restrictions on the permutations a (that they be fixed point free and involutions), 
Baik and Rains Q have shown that the limiting laws for the length of the longest increas- 
ing/decreasing subsequence are now the limiting distributions Fi and F4 for the scaled largest 



eigenvalue in GOE and GSE, see (4T) and (14. 3| ). Generalizing to signed permutations and 
colored permutations the present authors and Borodin [^, |ll| showed that the distribution 
functions of the length of the longest increasing subsequence involve the same F2. 

Johansson |23] showed that the shape fluctuations of a certain random growth model, again 
appropriately scaled, converges in distribution to F2. (This random growth model is intimately 



related to certain randomly growing Young diagrams.) In subsequent work, Johansson |24| 



showed that the fluctuations in certain random tiling problems (related to the Artie Circle 



Theorem) are again described by i<2. Finally, Johansson [p4| and the present authors |44| have 
considered analogous problems for random words and have discovered various random matrix 
theory connections. 
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